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ABSTRACT 

This  report  gives  details  of  a  five-component  strain-gauge  balance  that  has  been  developed  to 
measure  flow-induced  loadson  models  in  theDSTO  water  tunnel.  The  loads  are  very  small  and 
the  balance  was  designed  to  measure  side  and  normal  forces,  as  well  as  rolling,  pitching  and 
yawing  moments,  within  the  ranges  ±25  N ,  ±25  N ,  ±0.1  N  .m,  ±0.2  N  .m  and  ±0.2  N  .m  respectively. 
These  loads  are  at  least  2  to  3  orders  of  magnitude  smaller  than  those  typically  measured  on 
aircraft  models  in  the  low-speed  wind  tunnel  at  DSTO.  The  balance  has  not  been  designed  to 
measure  axial  forces.  Duetothesmall  loads,  it  was  necessary  to  use  semi -conductor  strain  gauges 
on  the  bal  ance.  The  f  i  ve-component  bal  ance  has  been  devel  oped  usi  ng  si  mi  I  ar  design  pri  nci  pi  es 
to  a  two-component  balance  developed  earlier  at  DSTO  for  use  in  the  water  tunnel. 
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Development  of  a  Five-Component  Strain-Gauge 
Balance  for  theDSTO  WaterTunnel 


Executive  Summary 

The  Defence  Science  and  Technology  Organisation  (DSTO)  water  tunnel  in  Air  Vehicles 
Division  has  been  used  extensively  over  the  years  to  carry  out  aerodynamic  research 
investigations.  Thetunnel  has  primarily  been  used  to  observe  detailed  flow  patterns  over 
models,  such  as  aircraft,  missiles,  ships  and  submarines.  In  recent  times,  the  ancillary 
equipment  used  with  thetunnel  has  been  developed  and  improved,  enabling  new  types  of 
measurements  to  betaken.  As  part  of  this  upgrade,  a  low-range  five-component  strain- 
gauge  balance  has  been  developed,  enabling  very  small  flow-induced  forces  and  moments 
on  models  to  be  measured.  This  followed  the  successful  implementation  of  a  two- 
component  balance  into  thetunnel.  Both  balances  have  been  built  using  similar  design 
methods. 

The  new  balance  can  be  used  to  measure  side  and  normal  forces,  as  well  as  rolling, 
pitching  and  yawing  moments,  within  the  ranges  ±25  N,  ±25  N,  ±0.1  N  .m,  ±0.2  N  .m  and 
±0.2  N  .m  respectively.  These  loads  are  at  least  2  to  3  orders  of  magnitude  smaller  than 
those  typically  measured  on  aircraft  models  in  the  low-speed  wind  tunnel  at  DSTO.  The 
balance  was  not  designed  to  measure  axial  forces,  since  such  forces  on  water-tunnel 
models  cannot  be  scaled  to  correspond  to  those  on  full -size  vehicles  in  air.  The  reason  for 
th  i  s  i  s  that  bou  nd  ary  I  ayers  on  w  ater-tu  n  nel  mod  el  s  are  I  ami  nar  w  hereas  those  on  f  u  1 1  -si  ze 
vehicles  are  generally  turbulent,  i.e.  there  are  different  drag  characteristics  for  the  two 
cases.  Semi-conductor  strain  gauges  have  been  used  on  the  new  balance,  as  for  the  two- 
component  balance.  The  balance  was  calibrated  manually  using  a  conventional  dead¬ 
weight  procedure,  whereby  different  types  of  forces  and  moments  were  systematically 
applied  to  the  balance  and  corresponding  output  voltages  from  the  five  channels  were 
measured,  enabling  the  calibration  relationships  to  be  determined.  For  a  first-order 
calibration,  the  standard  errors  were  found  to  be  less  than  0.1%  for  the  different  load 
components,  suggesting  that  loads  measured  by  the  water-tunnel  balance  will  be  of 
acceptable  accuracy.  Details  of  the  design,  manufacture  and  calibration  of  the  five- 
component  balance  are  given  in  this  report. 
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Notation 


[A  ] 
[E] 
e 


F X,  Fy,  F 


z 


f 

H 


Matrix  defined  by  equation  A 19. 

Matrix  defined  by  equation  A 18. 

Sumof  squares  of  residuals,  (used  when  calculating  calibration  coefficients  for 
a  strain-gauge  balance). 

Axial,  sideand  normal  forces  respectively  associated  with  the  balance  and  the 
model  coordinate  systems,  (N).  Positive  directions  are  given  in  Figure  3. 

N  umber  of  degrees  of  freedom  i  n  the  cal  i  bration  equations. 

Load  applied  to  a  strain-gauge  balance. 


H  Load  estimated  using  a  calibration  equation. 

K  Gauge  factor  of  a  strai  n  gauge,  K  =(AR/ R)/  (AL/ L) 

L ,  AL  Length  of  a  strain-gauge,  change  in  length  of  a  gauge  under  load  (m). 

Mx,My,Mz  Rolling,  pitching  and  yawing  moments  respectively  associated  with  the 
balance  and  the  model  coordinate  systems,  (N  .m).  Positivedirectionsaregiven 
on  Figure  3. 

N  Total  number  of  points  used  in  a  calibration. 

P  Index  of  summation. 


R,AR 

R 

se 

V  OUT,  V  IN 

x,  y,  z 


Electrical  resistance  of  a  strain-gauge,  change  in  resistance  of  a  gauge  under 
load  (Q). 

Voltage  ratio,  R  =I/0ut/I/in. 

Stand  ard  error,  used  to  assess  the  accu  racy  of  a  cal  i  brati  on . 

Output  and  input  voltages  respectively  for  a  channel  on  a  strain-gauge 
balance,  (V). 

Axes  for  the  balance  and  model  coordinate  systems.  Positive  directions  are 
given  in  Figure 3. 


Subscripts 

2,  3, 4, 5,  6  Refer  to  channels  2  to  6  on  the  five-component  strain-gauge  balance. 
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1.  Introduction 

Equipment  used  with  the  water  tunnel  at  the  Defence  Science  and  Technology  Organisation 
(DSTO)  has  recently  been  developed  and  improved  to  enablenew  types  of  measurements  to 
be  obtained.  As  part  of  the  upgrade,  a  technique  has  been  developed  to  measure  the  very 
small  flow-induced  pressures  on  the  surface  of  models  in  the  tunnel  (seeErm,  2000).  A  low- 
range  two-component  strain-gauge  balance  system  has  also  been  manufactured,  enabling  the 
very  small  flow-induced  normal  forces  and  pitching  moments  on  models  in  the  tunnel  to  be 
measured  (seeErm,  2006  a).  A  dynamic-testing  capability  has  also  been  developed,  enabling 
loads  to  be  measured  on  a  model  while  it  is  in  motion,  undergoing  a  specified  dynamic 
manoeuvre  (see  Erm,  2006  b). 

Following  the  successful  commissioning  of  the  two-component  balance,  a  five-component 
bal  ance  has  now  been  bu  i  1 1  usi  ng  si  mi  I  ar  desi  gn  methods.  O nee  agai  n,  semi  -cond uctor  strai  n 
gauges  have  been  used  on  the  balance.  The  new  balance  can  be  used  to  measure  side  and 
normal  forces,  as  well  as  rolling,  pitching  and  yawing  moments,  within  the  ranges  i25  N, 
i25  N ,  iO.l  N  .m,  i0.2  N  .m  and  i£>.2  N  .m  respectively.  The  balance  has  not  been  designed  to 
measureaxial  forces,  si  nee  scaled  tangential  skin-friction  drag  forces  on  water-tunnel  models 
are  not  representative  of  those  on  ful  I  -si  ze  ai  rcraft.  The  reason  for  thi  s  i  s  that  bou  ndary  I  ayers 
on  models  in  water  tunnels  aregenerally  laminar,  whereas  corresponding  boundary  layers  on 
full-size  vehicles  are  predominantly  turbulent.  By  restricting  the  balance  to  the  above  five 
components  (no  capability  to  measureaxial  forces),  it  isgenerally  notpossibleto  measure  lift 
forces  on  models,  which  are  normal  to  the  direction  of  thefree-stream  velocity.  Thisisamajor 
shortcoming  of  this  and  other  water-tunnel  balances.  However,  the  balance  significantly 
increases  the  usefulness  of  the  tunnel,  since  it  is  still  possible  to  measure  a  wide  range  of 
different  types  of  loads  on  models,  while  simultaneously  capturing  images  of  the  flow  over 
the  models,  enabling  the  loads  and  the  flow  patterns  to  be  correlated  directly.  Details  of  the 
design,  manufacture  and  calibration  of  the  five-component  balance  are  given  in  this  report. 


2.  DSTO  Water  Tunnel 

The  DSTO  water  tunnel  was  manufactured  by  Eidetics  International  Incorporated1  and  is 
designated  Model  1520.  The  tunnel,  shown  in  Figure  1,  has  a  horizontal -flow  test  section 
380  mm  wide,  510  mm  deep  and  1630  mm  long.  It  is  a  recirculating  closed-circuit  tunnel  and 
there  i  s  a  free  w ater  su rface  i  n  the  test  secti  on .  The  si  de  w al  I  s  and  f  I  oor  of  the  test  secti  on  are 
made  from  glass  to  facilitate  flow-visualisation  studies.  Thefree-stream  velocity  in  the  test 
section  can  be  varied  between  Oand  0.6m/  s.  The  contraction  upstream  of  the  test  secti  on  has 
an  inlet/  outlet  area  ratio  of  6:1.  There  are  six  dye  canisters  on  the  tunnel  that  can  be 
pressurised  with  air  to  force  dye  through  plastic  tubes  to  selected  locations  on  a  model  for 
flow-visualisation  studies.  There  is  a  suction  pump  on  the  tunnel  that  can  be  used  to  suck  or 
blow  water  through  models,  such  as  when  studying  the  effects  of  intake  flows  on  the 
aerodynamic  performance  of  an  aircraft.  Mod  els  are  mounted  on  a  C -strut  so  that  the  required 


1  Now  called  Rolling  Hills  Research  Corporation,  420  N.  Nash  St.,  El  Segundo,  CA  90245,  USA. 
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centre  of  rotation  of  a  model  is  at  the  centre  of  the  imaginary  circle  formed  by  the  strut, 
ensuring  that  all  angular  motion  of  the  model  is  about  this  point.  The  model  supporting 
system  isattached  tothetop  of  the  test  section  by  a  hinge,  so  that  a  model  can  belowered  into 
thetest  section,  or  removed  from  the  test  section,  as  required.  Further  detailsof  thetunnel  are 
given  by  Erm  (2006a). 


(b) 


Figure  1:  Eidetics  Model  1520  water  tunnel,  (a)  diagrammatic  view  of  the  tunnel, 
(b)  photograph  of  thetunnel 
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3.  Strain-Gauge  Balance 

A  sensitive  five-component  strain-gauge  balance,  shown  in  Figure  2,  has  been  designed, 
manufactured  and  calibrated  to  measure  side  and  normal  forces,  as  well  as  rolling,  pitching 
and  yawing  moments.  The  construction  of  the  balance  is  based  on  the  traditional  cantilever- 
beam  type  of  design,  and  consists  of  a  system  of  precisely-machined  flexure  members  upon 
which  strain  gauges  are  attached.  The  dimensions  of  the  flexure  members  and  the  likely 
deflection  of  the  balance  under  load  were  determined  using  a  finite-element  analysis.  The 
balanceisan  integral  unit,  machined  out  of  stain  I  ess  steel  rod,  having  a  specification  A I  SI  630 
Condition  H  480.  A  model  isattached  to  the  balance  at  its  tapered  end,  and  theother  end  of  the 
balance  is  mounted  on  a  tunnel  supporting  system,  such  as  a  roll-pitch-yaw  rig.  When  the 
model  is  loaded,  thesurface  strains  in  theflexure  members  cause  theattached  strain  gauges  to 
deform.  This  causes  the  resistances  of  the  gauges  to  change,  which  i  n  turn  affects  the  output 
voltages  of  the  gauges.  Loads  on  a  model  can  then  be  determined  using  balance  calibration 
relationships  (see  Section  6).  Basic  principles  of  balance  construction  and  operation  aregiven, 
for  example,  by  Edwards  (2000). 

The  flow-induced  loads  on  aircraft  models  in  the  water  tunnel  are  small  by  conventional 
wind-tunnel  standards,  and  the  balance  has  been  designed  to  measuresideand  normal  forces 
within  the  range  ±25  N ,  rolling  moments  within  the  range  ±0.1  N  .m,  and  pitching  and  yawing 
moments  withi  n  the  range  ±12  N  .m.  The  load  ranges  were  determi  ned  by  seal  i  ng  force  and 
moment  data  measured  on  full-sizecombat  aircraft  during  manoeuvres.  For  comparison,  the 
balances  used  in  theDSTO  low -speed  wind  tunnel  (LSWT)  are  typically  designed  to  measure 
normal  forces  and  pitching  momentswithin  the  ranges  ±3500  N  and  ±300  N  .m  respectively. 

The  relatively  low  free-stream  velocities  used  in  thewater  tunnel  wi  1 1  produce  corresponding 
smal  I  loads  on  the  models.  The  need  to  detect  such  smal  I  loads  necessitated  the  use  of  semi¬ 
conductor  strai  n  gauges,  havi  ng  a  resi  stance  of  1000  fi,  a  gauge  factor  of  145,  and  d  i  mensi  ons 
of  1.27  mm  by  0.15  mm  (active  area,  length  by  width).  The  gauge  factor,  K,  is  given  by 
(AR/ R)/  (AL/ L),  whereR  is  the  resi  stance  of  a  strain  gauge,  AR  isthechange  in  the  resistance 
of  the  gauge  under  a  load,!  isthelength  of  the  gauge,  and  AL  isthechangein  length  of  the 
gauge.  Thegaugeswere  manufactured  by  PSI-TRON IX  Incorporated2  and  the  model  number 
of  thegauges  is  P01-05-1000.  The  balance  contains  4  gauges  in  each  of  5  different  channels,  to 
measure  the  5  different  load  components,  so  that  there  are  20  gauges  on  the  balance.  The 
positioning  of  the  gauges  on  the  flexure  members  for  each  of  the  5  channels  is  as  shown  in 
Figure 3.  For  the  rolling- moment  channel,  gauges  are  positioned  at  45°  on  the  upper  and 
lower  flexure  members  (but  not  the  sideflexure  members).  The  four  gauges  for  each  channel 
have  been  positioned  on  flexure  members  and  wired  together  in  a  Wheatstone-bridge 
configuration  to  minimise  the  output  response  of  each  bridgedueto  other  load  components, 
i.e.  cross  coupling.  However,  in  practice,  there  is  a  small  amount  of  cross  coupling 
(see  Figure  7)  and  the  different  load  components  are  computed  using  the  outputs  from  all 
5  channels. 


2  PSI-TRONIX  Incorporated,  3950  South  “K”  Street,  Tulare,  CA,  93274,  USA. 
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Figure  2:  Five-component  balance,  showing  the  main  dimensions 
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Gauges: 

2.  Sidewards  force 

3.  Normal  force 

4.  Rolling  moment 

5.  Pitching  moment 

6.  Yawing  moment 


Figure  3:  Five-component  strain-gauge  balance,  showing  the  position  ofthegauges  for  thedifferent 
channels 


The  right-handed  orthogonal  balance  coordinate  system  isshown  in  Figure  3.  Thecoordinate 
system  remains  fixed  to  the  balance,  with  the  origin  located  on  the  longitudinal  axis  of  the 
bal  ance  at  the  geometr  i  c  centre  of  the  gauges.  The  d  i  recti  ons  of  the  correspond  i  ng  forces  and 
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moments  are  as  show n ,  w  i th  the  posi ti  ve  d i  recti ons  of  the  moments  determi  ned  by  the  right- 
hand  screw  rule  about  the  respective  axes. 

The  strain  gauges  have  been  glued  onto  the  balance  using  a  bonding  material  designated 
M-BON  D  6003,  which  isthesame  as  that  used  for  conventional  strain  gauges.  The  gauges  and 
the  connecting  leads  have  been  waterproofed  by  coating  them  with  a  compound  designated 
M -coat  C3,  a  solvent-thinned  (naphtha)  RTV  (room  temperature  vulcanising)  silicone  rubber. 
For  additional  waterproofing,  the  gauges  and  leads  were  also  covered  with  a  thin  sleeve  made 
from  si  I  i  cone  ru  bber.  1 1  was  i  mportant  to  use  a  thi  n  si  eeve,  si  nee  previ  ousl  y  it  was  found  that  a 
thick  sleeve  stiffened  the  balance  and  caused  it  to  exhibit  hysteretic  behaviour  during 
oscillating  loading.  Figure  4  shows  the  gauges  glued  onto  the  balance  prior  to  waterproofing. 
The  gauges  werefound  to  beslightly  sensitive  to  both  temperature  and  light,  butthisdid  not 
haveany  significant  effect  on  measured  data.  When  gauges  are  wired  in  a  Wheatstone-bridge 
configuration,  theeffects  of  temperature  on  thedifferent  gauges  ideally  nullify  each  other,  so 
that  there  is  no  resultant  change  in  bridge  output  voltage,  provided  that  all  gauges  are 
subjected  to  the  same  temperature.  The  temperature  of  the  water  in  the  tunnel  was  very 
stable,  varying  typically  by  less  that  0.2°C  throughout  a  day,  and  the  balance  was  mounted 
deep  inside  a  model  when  testing,  well  protected  from  strong  light  sources  used  when 
visualising  the  flow. 


cm 

Figure  4:  Five-component  strain-gauge  balance,  showing  gauges  and  leads,  prior  to  waterproofing 


4.  Signal-Conditioning  System 

An  existing  six-channel  signal -conditioning  system  was  modified  and  recommissioned  for  use 
with  the  water-tunnel  five-component  balance.  The  system  can  be  used  to  null,  amplify  and 
filter  output  voltages  from  the  Wheatstone  bridges  on  the  balance,  before  the  signals  are 
sampled  by  a  PC-based  data-acquisition  system  (see  Section  5).  A  simplified  circuit  diagram  of 
one  of  the  bridges,  the  power  supplies,  a  nulling  unit,  an  amplifier,  a  filter,  and  the  data- 
acquisition  system,  are  shown  diagrammatical  I y  in  Figure  5. 


3  M-BOND  600  and  M-Coat  C  are  manufactured  by  Vishay. 
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Bridge  on 

strain-gauge  balance 


Figure5:  Simplified  circuit  diagram  ofabridgeon  the  balance,  together  with  ancillary  components 


An  adjustable  power  supply  is  used  to  providea  +5.0  V  DC  excitation  voltage  to  the  strain- 
gauge  bridges  and  the  nul  I  i  ng  units,  and  another  power  supply  is  used  to  providea±15  V  DC 
excitation  voltage  to  the  amplifiers  and  filters.  Each  channel  contains  a  bridge  nulling  unit, 
with  both  coarse  and  fine  adjustments,  which  can  be  used  to  null  the  output  voltage  from  a 
Wheatstone  bridge,  i.e.  the  bridge  output  voltage  can  be  set  to  0  V.  Each  channel  also  contains 
an  amplifier,  whose  gain  can  be  set  to  100, 200, 400, 500 1000  or  2000,  as  well  as  a  third -order 
Butterworth  low-pass  filter  whose  3  dB  cut-off  frequency  can  be  set  to  1,  2,  5, 10,  20,  50, 100, 
200  or  500  Hz.  The  filter  can  also  be  switched  out,  in  which  case  the  cut-off  frequency  of  the 
filter  is  about  50  kHz,  which  is  the  inherent  bandwidth  of  an  amplifier. 

When  nulling  a  bridge,  the  voltage  at  the  output  of  a  filter  (points  E  and  F  on  Figure  5)  is 
monitored  rather  than  the  voltage  at  the  output  of  the  bridge  (points  A  and  C  on  Figure  5). 
Thevoltageattheoutputof  thefi  Iter  can  be  read  from  a  dial  on  a  console,  or  alternatively,  can 
be  read  by  connecting  a  voltmeter  across  the  terminals  labelled  OUTPUT  on  the  console.  A 
sampling  program,  named  NULL  BALANCE,  based  on  LabVIEW4  software,  has  also  been 
developed  so  that  the  voltage  can  be  read  on  a  monitor  using  a  graphical  user  interface. 


5.  PC  and  Data-Acquisition  Card 

The  PC  used  has  a  Windows™  XP  Professional  operating  system,  a  3.2GHz  central 
processing  unit,  a  220  GB  hard  disk  drive  and  1GB  of  RAM . 

The  data-acquisition  card  interfaced  with  the  PC  was  manufactured  by  National 
I  nstruments™  and  the  product  code  is  N I  6013.  The  card  features  16  channels  (8  differential) 
of  16-bit  analogue  input,  a  68  pin  connector  and  8  lines  of  digital  input/  output.  Analogue 


4  LabVIEW™,  a  graphical  programming  language,  is  a  product  of  National  Instruments. 
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input  voltages  varying  between  -5V  and  +5V  can  be  sampled.  For  the  16-bit  card,  the 
resolution  of  the  sampled  voltages  is  10.0/ 216,  i.e.  10.0/ 65536  =  0.000153 V/ LSB 
=  0.153  mV/  LSB  (LSB  denotes  least  significant  bit). 


6.  Calibration  of  the  Strain-Gauge  Balance 


6.1  Balance  Calibration  Equations 

The  for  mat  of  the  cal  i  brati  on  rel  ati  onsh  i  ps  d  epend  s  u  pon  w  h  i  ch  parti  cu  I  ar  bal  ance  cal  i  brati  on 
model  or  set  of  equations  is  chosen  to  represent  the  data.  There  are  a  number  of  different  sets 
of  equations  that  can  be  used,  but  the  one  currently  used  at  DSTO  for  the  wind-tunnel 
bal  ances  assu  mes  that  bal  ance  output  vol  tages  aref  u  ncti  ons  of  the  cal  i  brati  on  coeff  i  ci  ents  and 
the  applied  loads  (see  Lam  1989,  Leung  &  Link  1999,  Blandford  2004).  This  model  isalsoused 
for  the  water-tunnel  balances. 

Ideally,  the  output  voltage  from  say  the  normal -force  channel  isonly  affected  by  the  normal 
force  applied  to  the  bal  ance,  and  similarly  for  other  channels.  However,  in  practice,  balances 
are  not  ideal  and  there  may  be  interactions  between  the  different  channels,  so  that  the  output 
voltage  from  each  channel  generally  depends  on  all  of  the  load  components  applied  to  the 
balance. 

For  the  five-component  balance,  the  calibration  is  always  described  by  five  different 
equations,  but  the  number  of  terms  in  each  equation  can  vary  depending  on  whether  the 
equations  are  first-order,  second-order,  third-order  and  so  on.  A  first-order  calibration  is 
described  inthisreport.Thefirst-orderequationswill  also  be  used  for  actual  measurements  in 
the  water  tunnel. 


The  calibration  procedure  (see  Section  6.2)  would  have  been  complicated  significantly  by 
usi  ng  a  second-  or  thi  rd-order  cal  i  brati  on .  1 1  wou  I  d  have  been  necessary  to  apply  two  or  more 
load  components  simultaneously  during  calibration  (such  as  loads  applied  in  they  and  z 
directions).  For  a  well-designed  balance,  second- and  thi  rd-order  terms  aresmal  I .  Forthefirst- 
order  calibration  used,  the  standard  errors  for  the  bal  ance  were  found  to  be  small,  being  less 
than  0. 1%  for  the  d  i  fferent  I  oad  components  (see  Secti  on  6.5),  show  i  ng  that  such  a  cal  i  brati  on 
is  adequate  for  the  type  of  testi  ng  to  be  undertaken  i  n  the  water  tunnel . 

The  five  first-order  equations  are  given  by 


R  2  =  C  2  2 W  2  C  2,3 H  3  +  C  2  4H 

R  3  =  C  3  2  +  C  3  3 H  3  +  C  3  4 H 

R  4  =  C  4  2  ^  4  3 H  3  +  C  4  4 H 

R  5  =  C  5  2  C  5,3 W  3  +  C  5  4 H 

R  6  =  C  6  2  +  6,3^  3  +  6,4^ 


4  +C2j5H  5  +C2j6^  6 
4  +C35H  5  +c36h  6 

4  +  ^  4,5^  5  +  ^  4,6^  6 
4  +C55H  5  +C56H  6 

4  +  C6;5H  5  +  c  6  6h  6 


(1) 

(2) 

(3) 

(4) 

(5) 
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R  2,  R  3,  R  4,  R  5  and  R  6  are  the  vol  tage  rati  os  for  channel  s  2  to  6,  i  .e.  the  output  vol  tages  from  the 
channels  divided  by  the  corresponding  input  voltages  to  thechannels(R  =U0ut/  Vm)-TheH 
terms  are  the  corresponding  applied  loads  and  theC  terms  are  the  calibration  coefficients. 
N  ote  that  there  are  no  terms  for  axi  al  forces  (subscri  pt  1)  for  the  five-component  balance. 

The  five  first-order  equations,  expressed  in  matrixform,  are  given  by 


■«2" 

^  2,2 

^  2,3 

C  2,4 

^  2,5 

^  2,6 

1 

l\J 

_ 1 

A  3 

^  3,2 

^3,3 

^  3,4 

^  3,5 

^3,6 

h3 

r4 

= 

^  4,2 

^  4,3 

^  4,4 

^  4,5 

^  4,6 

h4 

r5 

^  5,2 

^  5,3 

^  5,4 

^  5,5 

^  5,6 

Hs 

«6 

f  6,2 

^  6,3 

C  6,4 

^  6,5 

^  6,6  _ 

h6 

Equation  6can  be  expressed  in  a  simplified  matrixform  as  follows: 

[R]=[C][H]  (7) 

6.2  Acquisition  of  Calibration  Data 

The  balance  was  calibrated  manually  using  the  experimental  setups  shown  in  Figure 6. 
Figure  6a  shows  thesetup  used  for  applying  sideforces,  normal  forces,  pitching  moments  and 
yawing  moments,  and  Figure  6b  shows  the  setup  used  for  applying  rolling  moments.  The 
balance  was  clamped  onto  a  calibration  supporting  rig  (not  shown)  so  that  its  centreline  was 
horizontal .  A  cal  i  bration  cage  was  bolted  onto  the  model  end  of  the  balance  at  the  tapered 
adaptor  (see  Figures  2  and  3).  The  rol  I  angle  of  the  balance  and  the  attached  cage  was  set  to 
either  0°,  90 e,  lEKPor  270 e,  as  appropriate.  The  only  contact  between  the  balance  and  the  cage 
was  at  the  tapered  adaptor.  For  thesetup  shown  in  Figure  6a,  a  pan  was  suspended  from  a 
knife  edge  along  the  top  of  the  cage.  Forces  and  moments  were  applied  to  the  balance  by 
placing  weights  in  the  pan  when  it  was  positioned  atdifferentx  locations  along  the  balance. 
For  the  setup  shown  in  Figure  6b,  the  pan  was  suspended  fromarod  attached,  inturn, to  both 
sides  of  the  cage,  and  forces  and  moments  were  applied  to  the  balance  by  placing  weights  in 
the  pan  when  it  was  positioned  at  differenty  locations  along  the  rod.  The  maximum  forces 
and  moments  applied  to  the  balance  during  a  calibration  were  chosen  to  correspond  to  the 
range  of  loads  I  i  kely  to  be  experienced  by  a  model  duri  ng  subsequent  tests. 

Detai  I  s  of  the  I  oad  i  ng  sched  u  I  e  correspond  i  ng  to  the  cu rrent  cal  i  brati on  are  given  i  n  Tabl e  1. 
Weights  were  progressively  added  to  the  pan,  and  then  progressively  removed,  for  the  pan 
positioned  at  different  positions  as  shown.  Altogether  there  were  14  different  loading/ 
unloading  sequences,  each  using  15different  loads  (which  indudesOg).  Prior  to  commencing 
each  loading  sequence,  when  there  were  no  weights  in  the  pan,  the  output  voltages  from  the 
five  channels  were  nulled,  i.e.  the  bridge  output  voltages  were  set  close  to  0.0  V.  The  method 
used  to  null  the  bridges  is  described  in  Section  4.  For  each  of  the  15  loads  used,  the  output 
voltages  from  the  five  channels  were  sampled  20  times  at  0.1s  intervals  and  each  set  of  20 
samples  was  then  averaged  to  obtain  mean  values.  Altogether,  210  calibration  points  were 
taken  for  each  channel  (15  loads  for  each  of  14  loading  sequences). 
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Figure 6:  Experimental  setup  used  when  calibrating  the  five-component  balance. 

(a)  setup  for  applying  side  and  normal  forces,  pitching  and  yawing  moments, 

(b)  setup  for  applying  rolling  moments. 


T ablel:  L  oading  schedule  used  when  calibrating  the  five-component  balance 


Channel 

Gain  of 

Amplifiers 

Filter  Settings 
(Hz) 

Coordinates  of  Load 

Application  (mm) 

Loading  Sequence 

(g) 

1  (Axial  Force,  Fx)  (not fitted) 

2  (Side  Force,  Fy) 

500 

1 

(1) x  =-10 

(2) x  =0 

(3)  x  =  +10 

Oto  140 to  Og  in 

20 g  increments 

3 (Normal  Force, Fz) 

500 

1 

(1) x  =-10 

(2) x  =0 

(3)  x  =  +10 

Oto  140 to  Og  in 

20 g  increments 

4 (Rolling  Moment,  M  x) 

100 

1 

NJ  M 

X  X 

II  II 

S5  $5 
>  > 

II  II 

Oto  35 to  0g  in 

5g  increments 

5 (Pitching  Moment,  M  y) 

500 

1 

(1) x  =-10 

(2) x  =0 

(3)  x  =  +10 

Oto  140 to  Og  in 

20 g  increments 

6  (Y aw  i  ng  M  oment,  M  z) 

500 

1 

(1) x=-10 

(2) x  =0 

(3)  x  =  +10 

Oto  140 to  Og  in 

20 g  increments 

A  sampling  program,  named  SCAN,  based  on  LabVIEW™,  was  developed  to  facilitate  the 
calibration  process.  For  each  of  the  210  sets  of  calibration  points,  the  user  had  to  specify  the 
number  of  samples,  time  between  samples,  outputfilename,  applied  load,  excitation  voltage, 
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temperature,  gain  settings  and  filter  settings.  The  program  was  then  run  and  the  output 
voltages  from  thedifferent  channels  weredisplayed  on  the  monitor  in  graphical  form  as  they 
were  being  sampled.  The  voltages  were  written  to  a  spreadsheet  file  along  with  other 
i  nformati  on  specif  i  ed  by  the  user.  I  nput  cond  iti  ons  were  updated  for  the  next  set  of  cal  i  brati  on 
points  and  the  procedure  was  repeated.  Fourteen  individual  datafiles  were  created  for  the 
comp  I  ete  cal  i  brati  on . 

6.3  Plots  of  Voltage  Rati  os  vs  Loads 

For  thecal  i  brati  on  data,  plots  of  V  out/  V  m  vs.  Fy,  M  z,  Fz,  M  y  and  M  x  areshown  in  Figure  7.  V  out 
isthe  output  voltage  from  a  Wheatstone  bridge  for  achannel  and  V|N  is  the  input  voltage  to 
the  bridgefor  that  channel  (\/in  is  nominally  5.0  V).  On  each  plot,  each  set  of  calibration  points 
for  the  different  channels  corresponds  to  two  loading/  unloading  sequences  (positive  and 
negative  values  of  Fy,  M  z,  Fz,  M  y  and  M  x).  Corresponding  points  for  the  loading  and  the 
unloading  phases  are  effectively  coincident  for  all  cases,  with  no  discernable  evidence  of 
hysteretic  behaviour.  The  voltage  ratios  vs  the  appl  ied  loads  of  al  I  5  channels  on  the  balance 
behave  linearly. 

6.4  Evaluation  of  Calibration  Coefficients 

The  calibration  coefficients  were  determined  from  the  discrete  appl  ied  loads  and  associated 
balance  output/  input  voltages  using  the  least-squares  regression  method  proposed  by 
Ramaswamy  et  al.  (1987).  Details  of  the  procedure  used  are  given  in  AppendixA. 

The  first-order  calibration  coefficients  are 


^  2,2 

^  2,3 

^  2,4 

^  2,5 

^  2,6 

'  0.22793 

-0.00962 

-0.38905 

0.36854 

0.16638' 

^  3,2 

^3,3 

C  3,4 

^  3,5 

^  3,6 

0.01056 

0.23426 

0.00390 

0.03656 

0.03167 

^  4,2 

^  4,3 

^  4,4 

^  4,5 

^  4,6 

= 

0.00635 

-  0.00248 

32.86634 

-  0.24769 

0.00160 

^  5,2 

C  5,3 

C  5,4 

^  5,5 

C  5,6 

-  0.00277 

0.00502 

-0.11615 

10.19147 

-0.34131 

f  6,2 

^6,3 

^  6,4 

^6,5 

C6,6_ 

-  0.00413 

0.00104 

0.09769 

0.53388 

10.09930 

6.5  Checking  Accuracy  of  Calibration 

When  calibrating  a  balance,  loads  are  appl  ied  to  the  balance  and  corresponding  output/  input 
vol  tage  rati  os  are  measu  red .  To  check  the  accu  racy  of  the  cal  i  brati  on,  the  reverse  proced  u  re  i  s 
used,  whereby  the  calibration  equations  are  used  to  estimate  loads  corresponding  to  the 
voltage  ratios  measured  during  the  calibration.  For  an  ideal  calibration,  the  applied  and 
estimated  loads  should  be  the  same,  butin  practice  small  discrepancies  do  occur.  Provided 
that  [C  ]  is  a  non-singular  square  matrix,  as  it  is  for  a  fi  rst-order  cal  i  brati  on,  then  [C  ] 1  (inverse 
of  [C  ])  can  be  formed  and  equation  7  can  be  rearranged  to  give 

[H]=mR]  (9) 

enabling  the  loads  to  be  estimated  directly. 
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-1.5  -1.0  -0.5  0.0  0.5  1.0  1.5 

Fy  (N) 


Figure  7a,  b,  c:  V  oltage  ratios  for  the  five  channels  vs  Fy  for  x  =  -10,  0,  +10  mm 
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-0.012  -0.008  -0.004  0.000  0.004  0.008  0.012 


Mz(N.m) 


Figure  Id,  e,  f:  V  oltage  ratios  for  the  five  channels  vsM  z  for  x  =  -10,  0,  +10  mm 
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-1.5  -1.0  -0.5  0.0  0.5  1.0  1.5 

Fz  (N) 


Figure  7g,  h,  i:  V  oltage  ratios  for  the  five  channels  vs  Fz  for  x  =  -10,  0,  +10  mm 
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-0.012  -0.008  -0.004  0.000  0.004  0.008  0.012 


/Wy(N.m) 


Figure  7j,  k,  I:  V  oltage  ratios  for  the  five  channels  vsM  y  for  x  =  -10,  0,  +10  mm 
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-0.015  -0.010  -0.005  0.000  0.005  0.010  0.015 

Mx(N.m) 

Figure  7m:  Voltage  ratios  for  the  five  chan  nets  vsM  xforx  =  +46.4,  y  =±40  mm 


A  statistical  indicator  commonly  used  to  assess  the  accuracy  of  the  calibration  isthe  standard 
error  (seeLam  1989  and  Leung  &  Link  1999).  The  parameter  indicates  the  good  ness  of  thefit 
of  the  cal  i  brati  on  equati  on  and  the  uncertai  nti  es  i  n  the  coeffi  ci  ents  determi  ned  by  the  I  east- 
squares  method.  The  standard  error  for  thesideforce  (for  example)  isgiven  by  the  foil  owing 
expression. 


se2  = 


(10) 


where  se2  is  the  standard  error  for  the  side  forces,  H2P  is  an  applied  si  deforce,  H2P  isthe 
corresponding  si  deforce  estimated  using  the  calibration  equation,  N  isthe  total  number  of 
points  used  in  the  calibration  for  each  channel  ( N  =210),  P  is  an  index  of  summation  and  f  is 
the  number  of  degrees  of  freedom  in  the  calibration  equations.  (The  number  of  degrees  of 
freedom  isequal  tothenumber  of  calibration  coefficients  per  component,  i.e.  f =5for  the  first- 
order  cal  i  brati  on  equations).  The  correspond  i  ng  termsfor  theother  load  components  (normal 
forces,  rolling,  pitching  and  yawing  moments)  aresimilarly  defined.  To  convert thestandard 
errors  to  dimensionless  parameters,  they  are  divided  by  the  appropriate  maximum  design 
load,  either  a  force  or  a  moment. 

Thestandard  errors  were  found  to  be  less  than  0.1%  for  the  different  load  components,  which 
is  comparable  to  the  standard  errors  of  0.19%  obtained  for  the  Collins  balance  used  in  the 
LSWT  -seeBlandford  (2004).  This  suggeststhat  loads  measured  by  thewater-tunnel  balance 
w  i  1 1  be  of  acceptabl  e  accu  racy. 
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6.6  Loads  M  easured  by  the  Balance 

When  carrying  out  tests  in  the  water  tunnel  with  the  balance,  loads  can  be  computed  from 
sampled  voltages  using  equation  9. 


7.  Checking  Waterproofing  of  Balance 

T ests  w  ere  u  ndertaken  to  see  w  hether  the  strai  n  gauges  and  the  el  ectri  cal  I  eads  on  the  bal  ance 
were  affected  by  water.  The  balance  was  submerged  in  water  for  a  period  of  10  hours,  and 
during  this  time  the  output  voltages  from  the  five  channels  were  measured  at  1  minute 
intervals.  The  voltages  showed  no  unusual  behaviour  over  the  10  hour  sampling  period, 
indicating  that  the  waterproofing  agent  on  the  balance  was  effective.  Submersion  of  the 
balance  in  water  had  no  noticeable  effect  on  subsequent  readings. 


8.  Concluding  Remarks 

This  report  gives  details  of  the  design,  manufacture  and  calibration  of  a  five-component 
strain-gauge-balance  used  to  measure  flow-induced  loads  on  models  in  the  DSTO  water 
tunnel.  The  loads  are  small  by  conventional  standards  and  thebalancecan  measuresideand 
normal  forces,  as  well  as  rolling,  pitching  and  yawing  moments,  within  the  ranges  ±25  N , 
±25  N ,  ±0.1  N  .m,  ±0.2  N  .m  and  ±0.2  N  .m  respectively.  The  balance  has  not  been  designed  to 
measure  axial  forces.  To  measure  the  small  loads,  it  was  necessary  to  use  semi-conductor 
strain  gauges  on  the  bal  ance.  The  calibration  of  thebalancewasfound  to  be  accurate  to  within 
±0.1%  for  thedifferent  load  components,  and  the  waterproofing  of  the  bal  ance  was  found  to 
be  effective. 
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Appendix  A:  Evaluation  of  Calibration  Coefficients 

In  Section  6.1,  it  was  indicated  that  the  first-order  calibration  relationships  for  the  five- 
component  balance  were 


R  2  =  c 2  2H  2  +  c  Z3H  3  +  C 2AH  4+CZ5H5+C 26H  6  (A  1) 

R3  =  C3>2H2+C3>3H3  +  C34H4+C35H5+C36H6  (A  2) 

R  4=C  4  2  H  2  +  C  43H  3  +  C  4 4H  4  +  C  4 5H  5  +  C  4  gH  g  (A  3) 

R  5  =  C  5.2H  2  +  C  5>3H  3  +  C  5  4H  4  +  C  5  5H  5  +  C  5>6H  6  (A 4) 

R  6  =  C  g  2H  2  +  C  6,3^  3  +  C  6  4/"/  4  +  C  6,5^  5  +  ^  6,6 ^  6  (A  5) 


R  2,  R  3,  R  4,  R  5  and  R  6  are  the  voltage  rati  os  for  channels  2  to  6,  i  .e.  the  output  voltages  from  the 
channels  divided  by  the  corresponding  input  voltages  to  the  channels.  TheH  terms  are  the 
corresponding  applied  loads  and  theC  terms  are  the  calibration  coefficients. 

The  coefficients  appearing  in  these  equations  can  be  determined  from  the  discrete  applied 
loads  and  associated  voltage  ratios  using  a  least-squares  regression  model.  Various  types  of 
regressi  on  methods  have  been  reported  i  n  the  I  iteratu  re,  but  the  method  used  for  cal  i  brati  ng 
the  water-tunnel  balances  is  that  developed  by  Ramaswamy  et  al.  (1987).  This  method  was 
also  used  when  calibrating  the  wind-tunnel  balances  in  Air  Vehicles  Division.  Lam  (1989) 
gives  mathematical  details  of  the  method,  and  only  an  outline  of  the  method  will  be  given 
here-also  see  Leung  &  Link  (1999)  and  Blandford  (2004).  The  method  will  only  bedescribed 
for  the  first-order  calibration  equations,  butthe  same  principles  apply  for  equations  of  higher 
orders. 


For  th  i  s  method ,  the  cal  i  brati  on  coeff  i  ci  ents  are  establ  i  shed  w  hen  the  su  m  of  the  squares  of  the 
differences  between  the  measured  strai  n-gauge output  voltage  ratios  and  those obtai  ned  from 
a  calibration  equation  is  a  minimum.  That  is, 


N  r  , 

e2  =  ^  r  2,2^  2,P  +  C  2,3^  3,P  +  ^  2,4^  4,P  +  ^  2,5 ^  5,P  +  ^  2,6^  6,P  _  ^  2,P  ] 
P=1 

e3  =  £[c  3,2^  2,P  +  C  3,3^  3,P  +  C  3AH  4p  +  C  3  5H  5  p  +  C  3  6H  6  p  R  3  p  ] 

P=1 

N  j.  , 

eA  =  2jf  4,2^  2,P  +  ^  4,3^  3,P  +  ^  4,4^  4,P  +  ^  4,5^  5,P  +  ^  4,6^  6,P  _  ^  4,P  J 
P=1 

e5  -  ^  E-  5,2^  2,P  +  C  5,3^  3,P  +  ^  5,4^  4,P  +  ^  5,5^  5,P  +  ^  5,6^  6,P  _  ^  5,P  ] 
P=1 

n  r 

e6  =  ^  6’2^  2’p  +  C  6’3^  3’p  +  C  6’4^  4,P  +  C  6,5^  5,P  +  ^  6,6^  6,P  _  R  6,P  J 


(A  6) 
(A  7) 
(A  8) 
(A  9) 
(A  10) 


are  al  I  mi  n  i  mu  m,  w  here  fe,  e3,  e4,  es  and  eg  are  su  ms  of  squ  ares  of  resi  d  ual  s  for  channel  s  2  to  6, 
P  is  an  index  of  summation  and  N  is  the  number  of  calibration  points  for  each  channel 
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(N  =210  -see  Section  6.2).  The  procedure  involves  partial  differentiating  with  respect  to 
each  of  i  ts  coeff i  ci  ents  and  equati  ng  the  resu  I  tant  expressi  ons  to  zero,  and  I  i  kew  i  se  for  63  to  65. 
This  results  in  a  set  of  five  simultaneous  equations  for  each  of  the  five  channels  as  follows5: 

Zf2,2H2+C„H3+C2i4H4+C2,5H5+C2,6H6-/!2]H2  =  0 
Zf2,2H2+C2,3H3+C2,4H4+C2.5H5+C2,6H6-«2]H3  =  0 
y[C23H2-C23H.+C24H4+C2=H5+C2,H,j-R2]H1  =  0  (All) 

Zf  22H  2  +C2,3H  3  +C2>4H  ,+Cz5H  ,+Cz,H 6 - R2]  H 5  =  0 
Zf2,2H2+C2,3H3+C2,4H4+C2.5H5+C2,6H6-«2]H6  =  0 

Zf3.2H2+Cs3H3+C34H4+C3.5H5+CMH6-R3]H2  =  0 
Z[CwH2+c3.3H3+C3,4H4+C3,5H5+Ci6H6-R3]H3  =  0 
IM  2+C33H3+C34H4+C35H5+C36H6-R3]H4  =  0  (A  12) 

Z[C3-2H2+C3,3H3+CS4H4+C3,5H5+C3.6H6-R3]H5  =  0 
Z[CwH2+c3.3H3+C3,4H4+C3,5H5+Ci6H6-R3]H6  =  0 

^  ]p  4,2^  2  +  C  4,3^  3  +  C  4,4^  4  +  ^  4,5 ^  5  +  ^  4,6 ^  6  —  ^  4]  ^  2  =  ^ 

[C  4,2^  2  +  C  4,3^  3  +  C  4,4^  4  +  ^  4,5^  5  +  ^  4,6^  6  —  ^4]  ^  3  =  ^ 

ZM  2  +  C 4  3H  3  +  C 4 4H  4  +  C 4 5H  5  +  C 4 qH  0  —  R4]h4  —  0  (A  13) 

^  [c  4,2^  2  +  C  4  3H  3  +  C  4  4H  4  +  C  4  5H  5  +  C  4  6H  6  —  R4]h5  =  0 

yi  [C  4,2^  2  +  C  4,3^  3  +  C  4,4^  4  +  ^  4,5^  5  +  ^  4,6^  6  —  ^4]  ^  6  =  ^ 

ZfwH2+CwH3+C5,4H4+c3,5H3+C3,6H6-R3]H2  =  0 

ZlC  5.2«  2  +' C 5,3H  3  +  C „H  4 + C 35H  5  +  C 36H  6  -  R 5]  H  3  =  0 

ZM  2+Ci3H3+Ci4H4+C55H5+C56H6-R3]H4  =  0  (A14) 

Z[C5.2H2+C3.3H3+C3,4H4+C3,5H5+C5,6H6-R5]h5  =  0 

Z[C5.2H2+C33H3+C5,4H4  +  C35H5+C5.6H6-«5]H6  =  0 


5  For  clarity,  the  subscript  P  on  R  and  H  is  omitted  in  the  remainder  of  this  report.  Whenever  the  summation 
symbol  is  used,  it  is  to  understood  that  the  sum  is  to  be  taken  over  the  range  P  =  1  to  N. 
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Z[C«H; 

>  +C 

6,3 H 

3+C 

6,+  ' 

^ + c 

6  ,sH 

5+C 

6,6^  6 

-««] 

El  2  - 

=  0 

?  +c 

6,3^ 

3+C 

6+  < 

^ + c 

6,5H 

5+C 

6,6^  6 

-fi6] 

\H  3~ 

=  0 

>  +  c 

6.3^ 

3+C 

6.4^  ^ 

i+C 

6.5^ 

5 

6,6^  6 

-««] 

h4- 

=  0 

l  +  C 

6,3  El 

3+C 

6,4+ 

^ + c 

6,5^ 

5+C 

6,6^  6 

-«6] 

\h5- 

=  0 

l  +  C 

6,3^ 

3+C 

6  4^  < 

^ + c 

6,5* 

5+C 

6,6^  6 

-k6] 

\h6- 

=  0 

(A  15) 


Equations  A 11  to  A 15  can  be  combi  ned  and  put  i  nto  matrix  notation  as  fol  lows: 

[e][c}t=[a] 

[C  ]T  is  the  transpose  of  the  matrix  of  cal  i  bration  coefficients  and  is  given  by 

re¬ 


fer  = 


2,2 

^  3,2 

E  4,2 

^  5,2 

^  6,2 

2,3 

^  3,3 

E  4,3 

^5,3 

^  6,3 

2,4 

^  3,4 

E  4,4 

^  5,4 

E  6,4 

2,5 

^  3,5 

+,5 

^5,5 

E  6,5 

2,6 

^  3,6 

^  4,6 

^  5,6 

E  6,6 

[E]  and  [A  ]  are  given  by 


[Eh 


ZH 

2  El  2 

ZH 

2h3 

ZH 

2H4 

ZH 

2  El  5 

ZH 

2h 

ZH 

3  El  2 

ZH 

3*3 

ZH 

3H  4 

ZH 

3H5 

ZH 

3H 

EH 

4  El  2 

ZH 

ro 

:c 

ZH 

4H4 

ZH 

4  El  5 

ZH 

4  H 

ZH 

5  El  2 

ZH 

5H3 

ZH 

5  El  4 

ZH 

5+ 

ZH 

5H 

ZH 

6^  2 

ZH 

6  El  3 

ZH 

6*4 

ZH 

6  El  5 

ZH 

eH 

(A  16) 


(A  17) 


(A  18) 


and 


[A]  = 


ZH  2R2 

ZH  2R3 

SH  2R4 

ZH  2R5 

ZH 

2r 

xh3r2 

ZH3R3 

sh3r4 

zh3r5 

ZH 

3r 

xh4r2 

sh4r3 

sh4r4 

sh4r5 

ZH 

4r 

zh5r2 

zh5r3 

sh5r4 

zh5r5 

ZH 

5R 

zh6r2 

zh6r3 

zh6r4 

zh6r5 

ZH 

6R 

(A  19) 


The  25  elements  comprising  [ E ]  and  the  25  elements  comprising  [A  ]  can  each  bedetermined 
from  the  known  applied  calibration  loads  and  the  corresponding  measured  voltage  rati  os,  so 
that  the  25  elements  in  [C]T  can  be  determined  by  making  use  of  equation  A 16.  It  is  not 
possibletodetermine[C]T  directly  from  thisequation,  and  theequation  must  be  rearranged  to 
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obtain  an  explicit  expression  for  [CF-  Provided  [£]  is  a  non-singular  square  matrix,  then  [E]1 
(inverse  of  [E])  can  be  formed  and  equation  A 16  can  be  rearranged  as  follows: 

[cV=mA]  (A  20) 


Once[C  F  is  known,  [C  ]  can  bedetermi  ned  by  si  mply  transposi  ng  [C  F  (equation  A 17)  to  give 


[C] 


c 

c 

c 

c 

c 


2,2 

^2,3 

^  2,4 

^2,5 

^2,6 

3,2 

^3,3 

^  3,4 

^3,5 

^  3,6 

4,2 

^  4,3 

^  4,4 

^  4,5 

^  4,6 

5,2 

^  5,3 

^  5,4 

^  5,5 

^5,6 

6,2 

^6,3 

^  6,4 

^6,5 

^  6,6 

(A  21) 
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